Microscopic simulation of superconductor-topological insulator proximity structures 
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We present microscopic, self-consistent calculations of the superconducting order parameter and 
pairing correlations near the interface of an s-wave superconductor and a three-dimensional topo- 
logical insulator with spin-orbit coupling. We discuss the suppression of the order parameter by the 
topological insulator and show that the equal-time pair correlation functions in the triplet channel, 
induced by spin-flip scattering at the interface, are of ± ipy symmetry. We verify that the spec- 
trum at sub-gap energies is well described by the Fu-Kane model. The sub-gap modes are viewed 
as interface states with spectral weight penetrating well into the superconductor. We extract the 
phenomenological parameters of the Fu-Kane model from microscopic calculations, and find they are 
strongly renormalized from the bulk material parameters. This is consistent with previous results 
of Stanescu et al for a lattice model using perturbation theory in the tunneling limit. 

PACS numbers: 73.20.-r, 75.70.Tj, 85.75.-d 



I. INTRODUCTION 

Fu and Kane showed that at the interface between 
a three-dimensional topological band insulator (TI) and 
an s-wave superconductor (S) forms a remarkable two- 
dimensional non-Abelian superconductor-i. It hosts Ma- 
jorana zero modes at vortex cores, as m. a Px + ipy 
superconductor—, but respects time-reversal symmetry. 
As argued in Ref. [l], the presence of superconductor 
induces a pairing interaction between the helical Dirac 
fermions at the surface of the topological insulator, and 
gaps out the surface spectrum. Then, the interface can 
be modeled elegantly by a simple matrix Hamiltonian in 
Nambu space (we follow the convention of Ref. ^ , 



iayAs 



(1) 



where k = {kx,ky) is the two-dimensional momentum 
in the interface plane, ai are the Pauli matrices, /is(k) 
is the surface Hamiltonian for the topological insulator 
describing the helical Dirac fermions^i^. 



/ls(k) = + Vs{(Jxky 



f^ykx)- 



(2) 



Fu and Kane also proposed to use S-TI proximity struc- 
tures to generate and manipulate Majorana fermions 
which obey non-Abelian statistics and are potentially 
useful for fault tolerant quantum computation^. This 
proposal and a few others that followed based on 
superconductor-semiconductor heterostructures^^— have 
revived the interest in superconducting proximity ef- 
fect involving insulating/semiconducting materials with 
spin-orbit coupling. More complex S-TI proxim- 
ity structures with ferromagnets"^^'^^ or unconventional 
superconductors^^ have been investigated. 

Given the importance of this proposal, it is desirable 
to understand to what extent the effective model Hpx 
holds, and what are the values of (As,/is,Us) for given 
materials. Answering these questions is crucial for future 
experiments designed to probe and manipulate Majorana 



fermions. As a first step in this direction, Stanescu et 
al considered a microscopic lattice model for the TI-S 
interface^. In this model, TI and S are described by a 
tight binding Hamiltonian defined on the diamond and 
hexagonal lattice respectively. The two materials are 
coupled by tunneling term in the Hamiltonian. These 
authors found that for small k, HfkO^) is valid but its 
parameters are significantly renormalized by the presence 
of the superconductor. This is supported by leading or- 
der perturbation theory in the weak coupling (tunneling) 
limit. They also discussed the induced p-wave correlation 
within the framework of perturbation theory. The p-wave 
correlation has also been noted in an analogous proxim- 
ity structure in two dimension between a quantum spin 
Hall insulator and a superconductor—. 

In this work, we consider S-TI proximity structures 
where S and TI are strongly coupled to each other, rather 
than being separated by a tunneling barrier. This is the 
desired, presumably the optimal, configuration to realize 
the Fu-Kane proposal, e.g. to achieve maximum value 
of As in HpK for given superconductor. In the strong 
coupling limit, the modification of superconductivity by 
the TI becomes important. This includes the suppression 
of the superconducting order parameter, the induction of 
triplet pair correlations by spin-active scattering at the 
interface, and the formation of interface states below the 
bulk superconducting gap. In order to accurately answer 
questions raised in the preceding paragraph for strongly 
coupled S-TI structures, one has to self-consistently de- 
termine the the spatial profile of the order parameter near 
the interface. 

Our work is also motivated by recent experimen- 
tal discovery that Copper-doped topological insula- 
tor Cua:Bi2Se3 becomes superconducting at a few 
Kelvin a^^'^^ . It seems possible then to combine such 
superconductors with topological insulator Bi2Se3 to 
achieve strong proximity coupling. Using Bi2Se3 and 
Cu2,Bi2Se3 as one of the examples, we set up micro- 
scopic, continuum models for the S-TI structures and 
solve the result Bogoliubov-de Gennes (BdG) equation 
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numerically. We first compute the superconducting or- 
der parameter as a function of the distance away from 
the interface. We then verify the validity of the Fu-Kane 
effective model and extract its parameters from the low 
energy sector of the energy spectrum. The emergence of 
Hfk will be viewed as the result of the "inverse prox- 
imity effect" , namely strong modification of supercon- 
ductivity by the presence of TI. This is in contrast to 
the previous viewpoint of pairing between surface Dirac 
fermions, which is a more proper description in the tun- 
neling limit. The spectral weight of these low energy 
modes (with energy below the bulk superconducting gap) 
are shown explicitly to peak near the interface but pen- 
etrate well into the superconductor. We will also show 
analytically that the induced triplet pair correlations are 
of Px ± ipy orbital symmetry, and systematically study 
their spatial and momentum dependence. Our results 
connect the phenomenological theory of Fu and Kanei 
to real materials. Our results for continuum models and 
strong coupling limit are also complementary to the re- 
sults of Stanescu et al^^- for lattice models and tunneling 
limit. 

In what follows, we first outline the formulation of the 
problem and then present the main results. Technical 
details on numerically solving the BdG equation are rel- 
egated to the appendix. 



II. MODEL AND BASIC EQUATIONS 

The band gaps of topological insulators are much larger 
than the superconducting gap of all weak coupling s-wave 
superconductors. For the purpose of studying the prox- 
imity effect between such superconductors and topolog- 
ical insulators, it is sufficient to describe the topological 
insulator using the low energy effective k • p Hamiltonian. 
Following Zhang et al^^, we model Bi2Se3 by 



i?T/(k) = 



/ M(k) Aik^ A2k^ 

M{k) A2k+ -Aik^ 

Aik, A2fc_ -M{k) 

V A2k+ -Aik,_ -M(k) 



Here k± = kx±iky, M{k) = M - Bikl~ B2{kl + kl), and 
/ is 4 X 4 unit matrix. The numerical values of the param- 
eters are obtained from first principle calculation o^^'^^ , 
M = 0.28 eV, Ai = 2.2 eVA, A2 = 4.1 eVA, Bi = 10 
eVA2, B2 = 56.6 eVA^. We work in basis {|1 t 
}, |1 I), |2 t), |2 ;)}, where 1 (2) labels the P1+ (P2+) 
orbitalii. Note that we have neglected the unimportant 
diagonal term eg (k) in Ref. |T3 which only slightly modi- 
fies the overall curvature of the band dispersion. We also 
keep the chemical potential /i as a tuning parameter. 

We consider a simple model of superconductor derived 
from a metallic state obtained by turning off the spin- 
orbit coupling {Ai = A2 = Q) in Hti and tuning the 
Fermi level well into the conduction band^-. The metal 
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FIG. 1: Schematic (not to scale) band diagrams in a 
superconductor-topological insulator (S-TI) proximity struc- 
ture. Ef is the Fermi energy of the metal described by Hm 
measured from the band crossing point. ^ is the chemical 
potential of TI measured from the band gap center. The su- 
perconducting gap is much smaller than the band gap of TI. 



Hamiltonian 

HM{k) = diag[M(k), M(k), -Af(k), -M(k)] - Efl, (4) 

with Ef > M . This mimics electron-doping the topolog- 
ical insulator^® or equivalently electrochemically shifting 
its chemical potential by applying a gate voltage^^. As 
shown in Fig. (TJ the valence band (band 1 with dispersion 
M(k) — Ef) is well below the Fermi level and remains in- 
ert as far as superconductivity is concerned. Next, within 
the framework of Bardeen-Cooper-Schrieffer theory, we 
assume attractive interaction between the electrons in 
the conduction band (band 2) near the Fermi surface de- 
scribed by the reduced Hamiltonian, 



H,^t = ^^^^(k)^^^(-k)A + h.c. 



(5) 



Here A is the superconducting order parameter, is 
the electron creation operator for orbital Z = 1,2 and 
spin a =t, 4-. The superconductor is then described by 

Hs=J2 ^li'^)HMik)i.M^iaik) + H,nt- (6) 

Note that Hs and Hti are in the same basis. 

This model is inspired partly by superconductor 
Cua;Bi2Se3. The transition temperature at optimal dop- 
ing X = 0.12 is Tc = 3.8K, which corresponds to a zero 
temperature superconducting gap A ~0.6meV— li^. The 
Fermi level is 0.25eV above the bottom of the conduction 
band, and the Fermi wave vector kf ^ 0.12A~^. The 
pairing symmetry of Cua;Bi2Se3 is to our best knowledge 
not clear at present. If it turns out to be a conventional s- 
wave superconductor, its mains features will be captured 
by Hs above with suitable choice oi Ef and A. 

Now consider a proximity structure consisting of a su- 
perconductor at z < d and a topological insulator at 
z > d (Fig. [T]). The interface at z = d is assumed to be 
specular, so the momentum k|| = {kx,ky) parallel to the 
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interface is conserved. The Hamiltonian for the whole 
system 

H = / (ik||d2;|^ vL(k||,z)[/io - /^(2:)]7/'l„(k||,z) 

- XI , z){hQ + ^(z)]v4„(k|| , z) 

+ A(z)V'^^(k|| , z)V'2i(-k|| , z) + h.c. 

+ Ai(z)[V'J^(-ia^)V'2t + ^u(*9^)V'2i + /i-c] 

+ yl2(z)[V'l^fc_ 1/^24 + V'ufc+V'2t + /i-c.]}. (7) 

Here /io(k||,92) — M — Bidl — -62^^, m(-z) and ^^(2;) are 
piece- wise constant, 

^l{z)^Ef0{d-z) + ^l6{z-d), (8) 
A,(z) = ^ = l,2 (9) 

in terms of the step function 9. The order parameter 
obeys the gap equation 



A(z)=5(z) / dk||(V'2t(k|h^)^2i(-k||,z)). 



(10) 



We assume g{z) = gO{d — z), the couphng constant g 
determines the bulk gap. 

To self-consistently solve Eq. (I?]) and (ITU|) . we intro- 
duce Bogoliubov transformation 



V'i<T(k||,z) = XM„,;^(k||,z)7„,k|| -t-<_,<,(k||,z)7;^ 



ra,k|| 



(11) 

to diagonalize H as 

n^Eg+ I dk|| X e„(fc|| 7„,k„ , (12) 

where £'g is the ground state energy, and 7^ ^.^^ is the cre- 
ation operator of Bogoliubov quasiparticles with energy 
e„(fc||). The wave function u and v satisfy the following 
Bogliubov-de Gennes (BdG) equation. 



iJB(k||,z)0„(k||,z) = e„(fc||)(/)„(k||,z). 
Here, the BdG Hamiltonian 



Hb 



f ho- d-cr 

d • <T —Hq — ^ — Ai(Tj^ 

^-/lo d-cr* 

\ A* iffy d- (T* /i + ho 



(13) 



(14) 



and the wave function (dropping the arguments) 

(t>n = iUn,lt,Un,ll,Un.2^,Un.2i,Vn,l'\-,Vn,li:^rL,2t,Vn.2i)'^ ■ 

(15) 

The vector d(k|| , z) is defined as 

= Ai{z)kx, dy = Ai{z)ky, d^ = A2{z){-idz). (16) 



Other quantities such as ho{k\\,z), fJ.{z), and A(z) are 
defined above. In terms of the wave functions, the zero 
temperature gap equation becomes 

Mz)=9iz) / dk|| XM„,2t(k||,2;X^2i(-k||,2), (17) 

where the summation denoted by prime is restricted to 
< en < ojD with UJD being the Debye frequency. 

We will exploit a particular symmetry of the BdG 
Hamiltonian to simplify calculations. Define the polar 
angle ifk for the in-plane wave vector k|| , 



kx ^ky ^ ^ 



(18) 



Then the BdG Hamiltonian for arbitrary {kx,ky) is re- 
lated to that for (fc^ — fc|| ,ky — 0) by unitary transfor- 
mation 

U\^)HB{kx,ky)U{k\\) = HBik\\,0). (19) 

Here ?7 is a block diagonal matrix, 

[/(k||) = diag[e-*'^^*,e-*'^-*,e"^^^,e"^-^]. (20) 

Thus, the eigen energy e„ only depends on the magnitude 
of k|| . Once the wave function for 95^ = is known, the 
wave function for ipk £ (0, 2tt) can be obtained by simple 
unitary transformation. 

We solve the matrix differential equation by con- 
serving it into an algebraic equation, following the treat- 
ment of superconductor-ferromagnet structure by Halter- 
man and Valla^S. The whole S-TI proximity structure is 
assumed to have finite dimension L in the z direction. 
The superconductor occupies the region < z < d, while 
the topological insulator occupies d < z < L. Hard wall 
boundary conditions are enforced at the end points, z — 
and z ~ L. The exact boundary conditions at the end 
points only affect the local physics there, provided that 
the boundaries are sufficiently far away from the S-TI 
interface. We expand the wave function and order pa- 
rameter in Fourier series-^. 



Un.la{z) 


'm 


(21) 


Vn,lcriz) 


^^V^nm^m{z), 


(22) 


A{z) 


X A„ (l)m{z), 

771 


(23) 


(t)in{z) 


= \^2/ Lsm{k^z). 


(24) 



The integer m ~ 1,2, ...,N labels the quantized longitu- 
dinal (along z) momentum km = mn/L. The cutoff N is 
chosen as^ 



(25) 



By expansion Eq. (|21l) - ([23|) . the BdG equation becomes 
an 8iV X 8iV matrix equation. With a reasonable guess 
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of the order parameter profile, the eigen energies and 
eigen wave functions are obtained by solving the matrix 
eigen value problem. Then a new order parameter profile 
is computed from the gap equation. The procedure is 
iterated until convergence is achieved. Relevant technical 
details can be found in the appendix. 

To analyze the spectrum of the system, it is convenient 
to define the retarded Green's function 

Gf,(k||,z,i) = -*0(O({Vi.(k||,z,t),7At^(k||,z,O)}) (26) 

where the time-dependent field operators are in Heisen- 
berg picture. For given k|| and the spectral functions 
are defined as 

Ni,m,z,uj) = -ImG«(k||,z,c^), (27) 
7V(k||,z,w) = ^7V,,(k||,z,c.). (28) 

la 

In terms of the wave functions and eigen energies, 

iV,„(k||, z, w > 0) = ^ lu„,i,(k||,z)p5(w - e„). (29) 

n 

We also introduce the equal-time pair correlation func- 
tions for the conduction electrons 

i^a/3(k||,2;) = (?A2a(k||,z)V'2/3(-k||,z)). (30) 

For example, at zero temperature we have 

^tt(k||,^) = Xl""'2t(k||,^X,2t(-k||,z), (31) 

n 
/ 

^u(k||,^) = Xl"«^2i(k||,zX^2i(-k||,z). (32) 

n 

Triplet components of F will be induced near the S-TI 
interface by spin-active scattering2^. 

III. THE ORDER PARAMETER 

First we present the spatial profile of the super- 
conducting order parameter ^{z) after the conver- 
gence is achieved. In all following calculations, Ef is 
fixed at 0.4eV, which is modeled after optimally doped 
Cua;Bi2Se3^. And the Debye frequency is set as ojd — 
O.lEf^. Fig. [5] shows an example with fj, = 0, L = 
BOOnm, d — 0.95L, and a bulk gap of 0.6meV as found 
in Cua;Bi2Se3. Going from the superconductor into the 
topological insulator, A first gets suppressed as the inter- 
face is approached before it drops to zero inside TI. The 
suppression is roughly 20% at the interface. Note that 
the fine wiggles of A in the simulation results are due 
to the finite momentum cutoff of the longitudinal mo- 
mentum kjn- As previously discussed by Stojkovic and 
Valls^^, the number of oscillations is ^ N/2, and the os- 
cillation amplitude vanishes in the bulk as N is increased. 
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FIG. 2: The superconducting order parameter A{z) near an S- 
TI interface at z = d = 0.95L. The superconductor occupies 
< z < d, and topological insulator occupies d < z < L. 
L — 300 nm, i-i—O, the bulk gap Ao =0.6meV. 



In this case, N is chosen to be 258 according to Eq. (|25p . 
So the matrix to be diagonalized is 2064 by 2064. 

Fig. [3]show the result for /i = 0, d = 0.9L, and a super- 
conductor with bulk gap Aq ~ 2.4meV. Since the coher- 
ence length is much smaller than the previous example, it 
is sufficient to consider L = 160nm, and correspondingly 
N — 138. The order parameter profile depends weakly on 
/I, as shown in Fig. |4]for a superconductor with bulk gap 
^ 5.2meV. From these examples, one observes that the 
length scale over which A is significantly suppressed does 
not scale with ^Oj the zero temperature coherence length 
of the superconductor. Rather it stays roughly the same, 
on the order of 30nm, as is varied over one decade 
from Fig. [2] to Fig. |4] (note the horizontal axis is z/L). 
This is not very surprising since is not the only length 
scale at play here. The interface represents a strong (as 
compared to Aq) perturbation that significantly distorts 
the bulk wave functions. The self-consistent microscopic 
BdG approach provides a reliable way to capture the de- 
tails of A(z) near the interface. 

It is illuminating to compare the proximity effect in 
S-TI structure with that in S-F structure^^, where F 
stands for a ferromagnetic insulator. The presence of 
F breaks time-reversal and spin rotation symmetry and 
significantly suppresses the order parameter. The sup- 
pression is sensitive to the spin mixing angle which is 
related to the band gap and exchange field of F— . In 
contrast, despite the spin-active scattering of electrons by 
TI which introduces spin-flips and spin-dependent phase 
shiftsS2j spin-orbit coupling is not pair breaking. The 
suppression of A near the interface is to a large extent 
due to the reorganization of local wave functions enforced 
by the boundary conditions at z = d for piece-wise po- 
tentials fJ-{z), Ai{z), g{z). It depends on for example how 
the wave functions decay inside the TI for given Ef and 
/i, and involves "high-energy" physics beyond the scale 
of A but below the scale of the band gap. To test this. 
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FIG. 3: The order parameter A{z) near an S-TI interface at 
z = d^ 0.9L. L = 160 nm, fi^O, Ao ~ 2.4meV. 
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FIG. 4: The order parameter profile for two different chemical 
potentials of the topological insulator, /i — — O.leV and /i = 
0.2eV. L = 160nm, Ao ~ 5.2meV. 



we have investigated the proximity effect between the 
same superconductor and a hypothetical ordinary insu- 
lator modeled by Hti with Ai = A2 = and the same 
band gap. The suppression of A by such an ordinary 
insulator turns out to be very similar. 



IV. THE INTERFACE MODE AND THE 
FU-KANE MODEL 

Next we analyze the energy spectrum of the system, 
£n{k\\), obtained from the BdG calculation. Take the 
case oi fj, = 0, L — 160nm, d = 0.9L, Ag ~ 5.2meV as an 
example. Fig. [5] shows the first several energy levels of 
the composite system versus the transverse momentum 
fc|| . There are many continuously dispersing modes at 
energies above the bulk gap. They are the usual Bogoli- 
ubov quasiparticles for different quantized longitudinal 



FIG. 5: The lowest few energy levels en(A:||). /i = 0, L = 
160nm, and the bulk superconducting gap Ao ~5.2meV. A 
well-defined interface mode is clearly visible at sub-gap en- 
ergies. Solid lines show a fit to the Fu-Kane model, with 
As = l.SmeV, Vs = 2.7eVA, and jj.s = 7.5meV. 



momenta. One also sees a series of avoided level cross- 
ings. At small fcji emerges a well-defined mode below Aq. 
We will identify it as the interface mode first discussed 
by Fu and Kane-. 

The Fu-Kane model Eq. ([T]) predicts the dispersion 



Eik) = v^A7F+(^a±M^ 



(33) 



We fit the very low energy portion of the spectrum to 
this prediction to extract the phenomenological parame- 
ters in the Fu-Kane model. The result is shown in Fig. 
[5l We find that, not surprisingly, A^ = l.SmeV which 
is much smaller than Aq — 5.2meV, and = 2.7eVA 
which deviates significantly from A2 = 4.2eVA predicted 
for the surface dispersion of TI. Moreover, /Zs = 7.5meV 
despite that the chemical potential of TI is /j, = 0. There- 
fore, our results show that the values of {As,Vs, ^J-s) are 
strongly renormalized by the presence of the supercon- 
ductor. This is consistent with the findings of Stanescu 
et al for weakly coupled S-TI structureaiS.. 

We have checked the validity of the Fu-Kane model for 
a variety of chemical potentials. Representative examples 
are plotted in Fig. [SI In each case, the sub-gap mode can 
be well accounted by the Fu-Kane model with suitable 
choice of parameters. While fig is always different from 
fj,, numerically we find it scales linearly with fj,. At the 
same time, A^ and show no strong dependence on fi 
for this set of parameters. To make sure that the sub-gap 
mode is indeed localized near the interface, we plot in Fig. 
[7] the z dependence of the spectral function N {k^\ , z , uj) . 
The spectral weight of the sub-gap mode is peaked near 
the interface and decays over a length scale ^ ^0 into the 
superconductor. Note that the spectral weight on the TI 
side (not shown in the figure) is finite, but it is much 
smaller in magnitude and decays very fast inside TI. 

We have carried out similar analysis for superconduc- 
tors with larger coherence length. Fig. [5] shows the evo- 
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FIG. 6: The dispersion of the lowest energy level for different 
fi (in eV). Other parameters are the same as in Fig. [5l L = 
160nm and Aq ~5.2meV. Fu-Kane model well describes the 
lowest energy mode. As /i is increased, As and Vs stay roughly 
the same, while Hs scales linearly with fj,. 




FIG. 7: The spectral function A'^(fe|| , z, cj) of the lowest energy 
level, oj = Emin, shown in Fig. [6] The interface is at 2: = 0.9L, 
Z/=160nm. The spectral function oscillates rapidly with z, so 
only its envelope is plotted. 

lution of the sub-gap mode with /i for Aq — 2.4meV. 
In this case, the values of {Ag^Vg, ii^) all varies with /i. 
Superconductors with larger ^0 E^nd smaller Aq are thus 
more sensitive to changes in fi and other microscopic de- 
tails near the interface. The exact values of the effective 
parameters in the Fu-Kane model in general depend on 
such microscopic details. 

V. TRIPLET PAIR CORRELATIONS 

It is well known that in heterostructures of s-wave 
superconductors, pairing correlations in other orbital 
channels, e.g. p-wave correlations, will be induced by 
scattering at the interface a^^'^^ . For example, inver- 



FIG. 8: The lowest energy level of an S-TI structure with 
L = 160nm, d — 0.91/, Aq — 2.4meV. /i is the chemical 
potential of the TI and measured in eV. 

sion/reflection symmetry (z ^ —2) is lost in an S- 
TI proximity structure, and the appearance of p-wave 
correlations seems natural from partial wave analysis. 
Moreover, scattering by a topological insulator is spin- 
active. The spin-orbit coupling inside a TI acts like 
a momentum-dependent magnetic field to flip the elec- 
tron spin and introduce different phase shifts for spin 
up and down electrons. The scattering matrix has been 
worked out by us previously---. Thus, a singlet s-wave 
Cooper pair can be converted into a pair of electrons 
in spin-triplet state at the S-TI interface. However, 
it is important to recall that by assumption attrac- 
tive interaction only exists (or is appreciable) in the s- 
wave channel. There is no binding force to sustain a 
triplet Cooper pair or a triplet superconducting order 
parameter. Similar (but different) pairing correlations in 
superconductor-ferromagnet hybrid structures have been 
extensively studied^^. The appearance of p-wave corre- 
lations in S-TI systems has been pointed out previously 
by Stanescu et al using a perturbative analysisi^. 

We focus on the equal-time pair correlation functions 
defined in Eq. pO|) . By exploiting the symmetry of the 
BdG Hamiltonian, Eq. (|19p . we are able to find ana- 
lytically the orbital structure of the triplet correlation 
functions. The unitary transformation Eq. ()20[) yields 

M2t(fc.,fcy) = ^i2t(fc|h0)e-^'^'=/2, 

U2i(fcx,fcj,) =M2i(fc||,0)e+^'^'=/2, 

«2i(fc.,fc^) =«2i(fc||,0)e-'''^'=/2. (34) 

Using these relations, we find 

^^tt(k||,z)-^^tt(fc|h^)e~*''N (35) 
n^(k||,z)=n^(fc||,z)e+^^'=. (36) 

Namely F^^ (Fn) hasp^-ipy (Px+iPy) orbital symmetry. 



7 




0.04 0.06 
k|| (1/A) 



0.6 




0.002 



0.004 0.006 
k||(l/A) 



0.008 



0.01 



FIG. 9: The imaginary part of triplet pair correlation function 
Ftt(fe||,2). The S-TI interface is at d = 0.9L. fj, = 0, L ^ 
160nm, Ao = 5.2meV. 



Finally, the remaining triplet correlation function 

(■02t(k|| , 2;)V'2i(-k|| , z) + V'24,(k|| , z)V'2t(-k|| , z)) (37) 

turns out to be zero. Note that the so-called odd- 
frequency paring correlations2i"i2^, which vanishes in the 
equal-time limit, are also interesting in S-TI structures, 
but we will not discuss their behaviors here. 

We find that F^^{k^\,z) is purely imaginary and iden- 
tical to Fii{k^^,z). The results ioi fi = 0, L = 160nm, 
d = 0.9L, Ao = 5.2meV are plotted in Fig. [H F-f^ 
vanishes at fc|| = as well as for large fc|| , namely when 

fc|| > -y/ {Ef + OJD + M)/i?2- This is consistent with lack 
of pairing in both limits. The behavior of F-f-^ for small fc|| 
is illustrated in Fig. [10]for ^ = 0, L = 300nm, d = 0.95L, 
Aq = 0.6meV. As comparison, we also plotted the singlet 
pair correlation function 

=Y1 "".2t(k|| , 2X,24,(-k|| , z) (38) 

n 

which is s-wave and purely real. 



FIG. 10; The imaginary part of _F-f^(fc|| , 2). fi — 0, L = 
300nm, d — 0.951/, Ao = 0.6meV. As comparison, the data 
points show the singlet pair correlation function F^i{k^^,z = 
0.9L)/3. 

low energy sector albeit with a set of renormalized pa- 
rameters. Our results are complementary to previous 
theoretical work on the proximity effectjiii^ and confirm 
the validity of the Fu-Kane model. 

We made a few simplifying assumptions in our calcu- 
lation. The superconductor is described by a two-band 
model with the valence band well below the Fermi level. 
Since only electrons near the Fermi surface are relevant 
for weak coupling superconductivity, we believe our main 
results are general. As idealizations, the chemical poten- 
tial, the spin-orbit coupling, and the attractive interac- 
tion are assumed to be step functions with a sudden jump 
at the interface. More elaborate and realistic models can 
be considered within the framework of BdG equations. 
For example, one can add a tunneling barrier between 
S and TI, or include a Rashba-type spin-orbit coupling 
term (due to the gradient of chemical potential) at the 
interface. We will not pursuit these generalizations here. 
Finally, the approach outlined here can be straightfor- 
wardly applied to study non-Abelian superconductivity 
in other superconductor-semiconductor heterostructures 
where spin-orbit coupling also plays a significant role^^— . 



VI. SUMMARY AND OUTLOOK 

In summary, we have investigated the proximity effect 
between an s-wave superconductor and a topological in- 
sulator using a microscopic continuum model. Strong 
coupling between the two materials renders the surface 
state of TI a less useful concept for this problem. Our 
focus has been on the various modifications to super- 
conductivity by the presence of TI. These include the 
suppression of the order parameter, the formation of in- 
terface modes below the bulk superconducting gap, and 
the induction of triplet pairing correlations. It is grati- 
fying to see the Fu-Kane effective model emerges in the 
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Appendix 

We follow the numerical scheme of Halterman and 
Vails to solve the matrix BdG equation^. The wave 
functions and the order parameter are expanded in the 
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orthonormal basis {(j)m{z)}, with m = 1,...,N. For ex- 
ample, function Un^i-[{z) is represented by TV numbers, 

/ It IT If It ^1 

Accordingly, each term in Hb is represented by a x iV 
matrix with the matrix elements given by 

/io(k||,a,) 5^^:{M - Bikl - B2kl) 

U{z) EfEjnm' + fJ-Fmm' 
A2{z)d, ^ A2G 

mm' 

Ai{z)k± A^k±F„un' 



These integrals can be evaluated analytically. Then the 
BdG equation becomes an 8A x 8A^ matrix equation. 
The gap equation can be rewritten as 



Am= g J dk|| ^ ^ Ji 



"nL'(k||X™"(-k||) 



2i 



m.m' ,m' 



n m' .7n" 



where 



Efnm' 
t^mm' 
Grnm' 
'Jm,m' ,m" 



(t>m{.z) (t>„i'{z)dz 



(t>m{z) 4>m'{z)dz 



(t)miz)d^(l)m'{z)dz 



(t>m{z) (t>ni'{z) (f)m"iz)dz 



The integral over k|| is first simplified to an integral over 
fc|l by the symmetry Eq. (|19l) and then evaluated numeri- 
cally with high momentum cutoff {Ep + cj/j + A/)/i?2- 
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